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ABSTRACT 

Using ’’Scenario Machine” we have carried out population synthesis of radio pulsar 
with black hole binaries (BH+Psr) in context of the most wide assumptions about 
star mass loss during evolution, binary stars mass ratio distribution, kick velocity and 
envelope mass lost during collapse. Our purpose is to display that under any suppo¬ 
sitional parameters of evolution scenario BH-l-Psr population have to be abundant in 
Galax y . It is shown that in t he all models inclu ding models evolved bv iHeger et alJ 
^l2002^ ■ IWooslev et alJ ll2002l) . lHeger et abl ll200,‘ll) expected number of the black holes 
paired with radio pulsars is sufficient enough to discover such systems within the next 
few years. 

Key words: stars: abundances - binaries: close - binaries: general - X-rays: binaries 
- pulsars: general - black hole physics. 


1 INTRODUCTION 

Binary radio pulsars with black holes (BH+Psr) reve¬ 
lation would be of fundamental significance for evidence of 
black holes exist ence and for investigation precision genera l 
relativity effects (iNaravan et al.iri99ll: iLinunov et al.lll994l . 
In such systems parameters of black holes - such as mass, 
Kerr parameter - will be measured with precision orders 
of magnitudes higher than indirect estimations of candi- 
dates in black holes masses in X-ray/BH-candidates binaries 
iBIandford fc Teukolsk\ll97RHBrumberg et al.ll97,'^^ . More¬ 
over, if mutual disposition is apt, it might be able to observe 
propagation of radio emission arbitrarily near to an event 
horizon. First accounts of possible number of BH+Psr bi¬ 
naries conducted ten years ago displayed the systems might 
be observed by the m odern radio-astronomical instruments 
iLinunov et ai jll994l . 

However total observed radio pulsars number have in¬ 
creased twice for the last 10 years and ran up to ss 1500 
but n one of them paired with blac k hole ||Manchester_et_alJ 

e^iTl20(^ 

mdowsk^^hl200^TAT^^ms^!atl200^ . Besides dur¬ 

ing these 10 years our conceptions of evolution of stars 
which are able to produce black holes appreciably change. 
In particular, considerations in favour of greater mass loss 
for these stars were obtained, detailed numerical computa- 
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tions considering new factors appeare d iHeger et al.ll2od^ : 
IWooslev et alJlioO^ : iHeger et alJl200,lll . 

We have carried out population synthesis of binary stars 
using ’’Scenario Machin e”. Description of it’ s working prin¬ 
ciples may be found in llLipunov et ahlllQQfil) . 

Since binary radio pulsars with black holes have to be 
generated by massive binary stars we relied on the observ¬ 
able statistics of the candidates for black holes paired with 
OB-stars (BH+OIU; _ 

In the paper llLipunov et alJll994l . it is supposed that 
any black hole paired with OB-supergiant must reveal it- 
self as a system Cyg X- 1 type. But, as it was shown by 
iKarpov fc LipunovI (1200 ih . powerful X-ray radiation is able 
to originate only if accretion disk around black hole has 
formed. Accordingly, in this paper, we term the Cyg X-1 
type systems as the subclass of such BH+OB binaries which 
has accretion disks. 

Note that Cyg X-1 type systems usually are not pro¬ 
genitors of BH+Psr binaries. More than 90% of them be¬ 
come BH+BH or single BH after merging of components 
during our calculations. The Cyg X-1 system has apprecia¬ 
ble chance of merging in the next stage - common envelope 
(CE) stage JBethe fc Brownl[l999ll . Ordinarily, BH+Psr bi¬ 
naries are originating from more wide systems which do not 
merge during CE stage or come through it. At that during 
BH+OB stage accretion disk does not form and observer is 
not able to see bright X-ray source. 

Evolutionary scenario which results in BH+Psr forming 
may be rough outlined as the next. We start calculation of 
massive binary system. When primary (more massive) star 
fills it’s Roche lobe mass transfer takes place and helium 





































2 V.M. Lipunov, A.I. Bogomazov, M.K. Abubekerov 


star remains instead of the primary star. As a rule black hole 
forms first, it’s companion - OB-star - is sufficiently sepa¬ 
rated from BH. At that in wide systems which can survive 
even after second mass transfer disk does not form (stellar 
wind velocity is too high near black hole and the moment of 
rotation of grasped matter is too little to form a disk). When 
second component fills it’s Roche lobe CE-stage begins and 
only wide systems (which number is greater than number of 
close binaries) produce BH-|-Psr (of course, anisotropic kick 
plays important role in this process, it leads not only dis¬ 
ruption of a system, but, if kick velocity and direction are 
apt, it can bound a binary during supernova explosion). In 
scenarios with high mass loss rate CE stage might not to be¬ 
gin, because mass of optical star is not enough greater than 
BH mass. Another scenario of BH-|-Psr creation is possible: 
when neutron star forms before black hole. At that time pul¬ 
sar has long and still purely explained history. Fortunately, 
part of these evolutionary tracks is small (see below). 

Note that evolutionary scenario contains quite a num¬ 
ber of key parameters poorly explained by theory (stellar 
wind magnitude, initial mass of a star which is able to form 
black hole, initial binary stars mass ratio distribution, kick 
velocity during relativistic stars formation, common enve¬ 
lope stage efficiency, part of mass falling into BH during 
collapse). Although i t is po ssible to reduce them significantly 
iLipunov et alll994Lll99^ in present work we suppose them 
as independent parameters. Qualitative influences of these 
parameters on stellar evolution scenario h ave been inves¬ 
tigated more than once in previous works Jde Jageilll980l: 
IShore et alJll994l: iLimmov et al.lll9^. 1 1 99711 Woosleve^^ 

l2002^ and it is possible to briefly delineate them in the fol¬ 
lowing way. 

Stellar wind magnitude essentially influences scenario 
for the two reasons. The hrst consists in that that spherically 
symmetric wind leads increasing of component separation. 
The second is in that total mass loss of a star by wind may 
cause change of it’s remnant type (it may produce neutron 
star instead black hole). 

Increasing initial mass Mmin of a star which is able to 
form black hole one decreases possible number of black holes 
due to Salpeter power law 0 . 

Initial binary stars mass ratio distribution is important 
because BH+Psr binaries have massive progenitors, so in 
case of more flat (a^— > 0) distribution 0 probability of 
BH+Psr systems birth decreases. 

Anisotropic kick during compact star formation is de¬ 
tailed investigated bv iLipunov et af] dlQQTll . Increase of kick 
velocity leads to total decay of a binary stars including rela¬ 
tivistic companion. It was shown that average kick velocity 
< V >« 150— 180 km s“^ accords to the observable number 
of binary radio pulsars. For the last few years two component 
kick velocity distribution with charact eristic velocities 90 
km s ~^ and 500 km s“^ were suggested i Arzoumanian et alJ 
1200 ( 11 . 

Effectiveness of CE stage is described by parameter 
CX-CE “ AEorfe, where AE^ = Egrav Ethevmai is 

the binding energy of the ejected envelope matter and 
AE„rh is the drop in systems orbital energy during spiral-in 
iShore et alJll99^ . Our population synthesis results shows 
very weak dependence on CE stage efficiency due to flat ini¬ 
tial distribution of major semi-axis of binaries (in logarith¬ 
mic scale). If one is decreasing acE close binary systems are 


originating from more and more wide binaries, which num- 
ber is not changin g on account of flat initi al distribution 
i jSlmre et alJIlQQ^ . We suggest acE ~ 0.5 iLiounov et alJ 
llOOill . that is satisfactory very close binary systems to form 
even in case of a very high stellar wind. Henceforward we do 
not investigate any dependencies on this parameter. 

Part of mass falling into BH during collapse is de¬ 
noted as kbh = Mbh/MpreSN, where Mbh is black hole mass, 
MpreSN - pre-supernova mass. It is important value because 
binary can experience possible decay concerned with quan¬ 
tity of mass ejected by a star during supernova explosion, i.e. 
the smaller kbh the smaller chance a system has to survive 
in the cataclysm. 

We emphasize that purpose of this work is not to hnd 
optimal model(s) of stellar evolution (another our work will 
be concentrated on it), but it is to display that radio pul¬ 
sar with black hole systems number in Galaxy have to be 
sufficient enough under any suppositional parameters of evo¬ 
lution scenarios to observe them within the next few years. 
So we have obtained multiple modelling using various ap¬ 
propriate parameters. 

At that, as in the work iLipunov et al.lli994l . we were 
calculating ratio of binary radio pulsars with black holes 
number to the total number of pulsars (practically the last 
number is the number of single pulsars). It allows to be 
saved from many selection effects concerned with our lack 
of knowledge of average polar pattern of the pulsars, their 
lifetime, magnetic field decay time, distribution of the ve¬ 
locity of the pulsars and their spatial distribution emerging 
at calculations of absolute number as long as we suppose 
that physical parameters of radio pulsars paired with black 
holes have no differences relative to average characteristics of 
the pulsars of the held. This natural suggestion is righteous 
because all of them have the same progenitors - massive 
(M > lOM©) stars and, as it was shown by our calculations, 
part of non-typical evolutionary tracks (when binary radio 
pulsars parameters might be dissimilar to parameters of sin¬ 
gle radio pulsars, for instance, due to effect of recycling) in 
BH+Psr binaries and any other radio pulsars (binary and, of 
course, single) is negligible (in most of models this quantity 
is smaller than 5% and never higher than 35%). 

Since BH+Psr systems have not been observed yet we 
suggest that observational limit is where 

1500 is rou nded number of obse rved radio pulsars - binary 
and single lATNF psr cat 


2 MODELS DESCRIPTION 

Common param eters for all of the models are 
iLipunov et alJTl99(1) : initial mass distribution (Salpeter 
function) 0, where Mi - initial mass of more massive star, 
mass ratio function 0, where q = M 2 /Mi - initial compo¬ 
nent’s mass ratio and a takes on a values 0 and 2, semi-major 
axis distribution 0 , where a is semi-major axis wit hin the 
limits of IORq < a < 10®E© iKraicheva al]ll981^ : 

f(M)=Mi-^-^^, (1) 


/(q) = 


( 2 ) 
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/(a) = i, (3) 

Kick velocity along with mass loss rate is one of the cru¬ 
cial and bad hxed parameters affecting population synthesis 
results. Higher kick leads to reduction of nu mber of binary 

S ns containing relativistic companions JLinunov et alJ 
. Although information about kick for neutron stars 
may be r eceived from the observable binary radio pulsars 
statistics llLvne fc Lorimeilll99'4 iHansen fc Phinne\^ll99'^ 
and from single p ulsars velocities, neu t ron stars kick is still 
under discussion llWiller ns et alJl2004 iMurphv et alJl20o3 
IPodsiadlowski etaHl^oi^ . 

In this work we have assumed the Maxwellian distri¬ 
bution of the kick velocity v during neutron star and black 
hole formation: 


(4) 

where no is characteristic kick velocity. 

Even less it is known of probable black hole kick, that 
is why we vary this quantity within the bounds of 0 
and 1000 km undoubtedly exceeding observable uncer¬ 
tainties, characteristic kick velocity absolute value for black 
holes also depends on part of mass lost during black hole 
formation: 


bh 


Vq = Vq 


MpreSN — Mhh 
Mhh 


( 5 ) 


where Mpresn - star mass before collapse, Mbh - black hole 
mass. 

The most important characteristic - mass loss rate of 
optical stars during evolution we have depicted by dint of 
five models - A, B, C, Wc and Wb. 

Scenario A has a weak stellar wind. Mass lo ss rate M 
during main sequence (MS) stage llde ,Tagerlll98ffl is: 


Mr^L/V^, (6) 

here L - star luminosity, Voo - wind velocity at infinity. 

For gia nts we take m aximum between and result 
obtained byE amer 3 (1198 j) : 

M H , (7) 

where R is stellar radius, M - it’s mass. 

And for red superg iants we take maximum bet ween 
and Reimers’s formula dKudritzki fc R,eimersllT984^ : 


M ~ LR/M, ( 8 ) 

Mass change AM in wind type A during one stage is 
no more than 0.1(M — Mcore), where M is star’s mass in 
the beginning of a stage, Mcore - it’s core mass. Mass loss 
during Wolf-Rayet (WR) star stage is parametrized as 0.1 • 
Mwr, where Mwr - maximum star mass during WR stage. 
We used fo r calculations of stellar wind type A core masses 
obtain ed bvIVarshavskii fc TutukovI dl97Rfl ; llben fc Tutuko\j 
lll98Slll98^. 


Scenario B uses c alculations of single-star evolution 
bv ISchaller et all (^O^). According to these calculations, a 
massive star loses most of its mass because of the action of 
stellar wind, down to ss 8 —IOMq before collapse, practically 
independent of it’s initial mass. 

In scenario C stellar e volution model is based on the re¬ 
sults of l Vanb evereni dlQQldl , which reproduce most accurately 
the observed galactic WR star distributions and stellar wind 
mass loss in massive stars. Mass loss by a star calculations 
were conducted if we used the next formula: 


AM = (M - Mcore), (9) 


where Mcore is stellar core mass dTHb -TTHH. If maximum 
star’s mass (usually it is initial mass of a star, but mass 
transfer in binary system is able to increase it’s mass over 
initial value) Mmax > ISM© mass of core in main sequence 
stage is determined as dTITbl. in giant and in supergiant 
stages as (llOt H. In Wolf-Rayet star stage, if Mwr < 2.5Mq 
and Mmax ^ 20Mq it is described as llOt Yl. if Mwr ^ 
2.5Mq and Mmax ^ 2 OM 0 as l[^), if Mmax > 20Mq - 
as IllOt 'l. 


^^core 


imM^ll, (a) 

— 3.051+4,21 Ig PI max —0.93(lg Plmax^"^ (/?) 

■ 0.83Myi, ’ ( 7 ) 

1.3 + 0.65(Mw'fl-2.4), (5) 

. S.OSM^lf, (e) 


( 10 ) 


Scenario C has high mass loss during the WR stage, it 
may reach 50% of a star mass or more here. Mass loss in 
other stages (MS, giant, supergiant) for stars with masses 
higher than ISM© (for less massive stars this scenario equals 
to type A wind) may reach « 30% of mass of a star. Total 
mass loss AM during all stages always is smaller than in 
scenario B and greater than in scenario A. 

W model is made in two types: with moderate (Wc) and 
with high (Wb) stellar wind. Pre-supernova, helium core and 
compact remnant masses according to the i nitial star mass 
to cal culate p arameter kbh were taken from lWooslev et alJ 
J20^, fig. lO. lHeger et al.l (l2002h . fig. 2. As it should be we 
calculate kbh as the next ratio: Mbh/MpreSN, where MpreSN 
is mass of pre-supernova star which has produced black hole 
with mass Mbh- We made our calculations in assumption 
that model Wc has type C wind, for model W b we used type 
B win d. This suggestion is quite right because lSchaller et alJ 
lll992^ made his calculation s using stellar wind obtained by 
iNieuwen huiizen fc de Jagej il99g) (thi s type of w ind was 
used by Iwoosle^^^^ |j20^ 1. jHe ger et alJ ^2^^2^ as high 
mass l oss type), IVanberorCT W9f ) - including mass loss 
rate by 
used by 


Hmnam^^j£oesterte I W9S ) (this type o f wind was 
^foos^n^^m^02^ . Heger et alJ l^2002^ as reduced 


mass loss type). 

Parameter kbh in A, B and C model is an adjusted con¬ 
stant value for all supernova explosions, in Wc and Wb mod¬ 
els it is various quantity dependent on initial mass of a star 
which is able to produce black hole. 

Finally it is necessary to introduce minimal mass of 
black hole Mmin (i.e. minimal pre-supernova star mass 
which is able to form black hole multiplied by kbh) which 
have been used for calculations in models A, B and 
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C. We vary this parameter within very wide bounds: 

2.5M0 ^ Mmin ^ lOM©. 


3 OBSERVATIONAL FOUNDATION 

To estimate probable number of binary pulsars with 
black holes there is a need to be guided by observable quan¬ 
tity of black holes candidates in our galaxy. From this view¬ 
point the nearest black hole paired with radio pulsar relation 
is Cyg X-1 type system - black hole paired with blue super¬ 
giant. Namely these binaries (BH-I-OB) are radio pulsars 
paired with black holes progenitors. As it is well known we 
observe only the one such source in the Galaxy for the time 
present - Cyg X-1. Although it i s not possible to speak about 
a stati stics, we suggest following ivan Paradiis fc McClintockI 
that this object is not a statistical ejection and to¬ 
tal number of such systems in Milky Way may reach a few. 
Similar candidates existence in the Large Magellanic Cloud 
- LMC X-1, LMC X-3 - convinces us of this. 

As we have already marked Cyg X-1 is not simply bi¬ 
nary consisting of black hole and blue supergiant but it is a 
very close binary system in which the regime of disk accre¬ 
tion on the black hole has been realized. It is not by chance - 
very low stellar win d velocity is necessary to form accretion 
disk ^LiDunovlll992^ : 

V < u,. « 320{Arj)imiT~^R~^{l + ta.n py^, (11) 

where rj - averaged over z-coordinate dynamic viscous coeffi¬ 
cient, m = Mx/Mq, Mx - relativistic star mass, Tio = T/10, 
T - orbital period in days, Rg = i?min/10®cm., Rmin - min¬ 
imal distance from the compact object up to which free Ke- 
plerian motion is still possible, /3 - accretion axis incline 
angle with respect to the radial direction. For black holes 
Rmin = 3Rg, where Rg = 2GMt,h/c^ ■ 

Inasmuch as stellar wind velocity rises with increasing 
of distance R fr om the normal star approximately as follow 
ide ,Iagerlll98rfl : 

V = V^(l-R/ay^^, ( 12 ) 

where a is characteristic radius of the star, accretion disk 
does not form in the most cases. Actu ally, velocities o f stel¬ 
lar winds are not precisely measured Jde ,Iageilll98fil . But 
GUl is quite good approximation for this work. A spheri¬ 
cally sy mmetric accretion could n ot set a bright source on 
the sky llKarnoy fc LiDunoyll200ih . Thus in the present ar¬ 
ticle we eyerywhere assume that the system Cyg X-1 type 
is a very close binary including blue supergiant with mass 
higher than 10 masses of the Sun in which disk accretion 
regime has been realized by the data iim . 

Let us estima t e l| llll and 11211 for the Cyg X-1. 
IPavis fc Hartmann! dlOS'lli obtained Uoo = 2300 ± 400 km 
s“^ for this system. It gives lower limit for wind velocity 
near black hole V ~ 1240 km s“^ I12II . Parameters in HII 
and GUl for Cyg X-1 has the next values (lAbubekerov et ^ 
l2004^ - R/a « 0.57, Tio = 0.56, m = 10, Rg = 0.1, 
tan/I « 0.1, also we suppose that (d?;)* « 1. So, critical 
velocity Vcrit ~ 1170 km s~^. Taking into account that ob¬ 
servational data have great uncertainties about Voo and o 
has approximate character (U < Vcr ~ ...) this coincidence 


can be considered as quite good. Hence our theoretical es¬ 
timation llTB and iHUl is in approximate agreement with 
observational data about accretion disk forming conditions 
in the Cyg X-1 system. 

We especially note that the Cyg X-1 is not the direct 
forerunner of the binary radio pulsar with black hole, most 
likely this system will merge after Roche lob e infill and com¬ 
mon envelope stage (iBethe fc Brownlll999ti . However, it is 
clear that Cyg X-1 type systems are very similar to BH-|-Psr 
progenitors which becomes apparent in correlation between 
their abundances during population synthesis (see below). 

During population synthesis we were picking out 
only the pulsars with black holes hav ing observable 
orbit al periods (less than 10 y ears) lLanib_£^X!amb| 
[]i)7(^__ |Majiclmsterfc Tavloil Il977t I.IohnstotT et al .1 Il992l: 
[ATNFDS^!ani200^i ~ it is may be hard to find larger peri¬ 
ods and the most known binary radio pulsar orbital period 
is less than 4 years. 


4 RESULTS 

The radio pulsar with black hole binaries quantity de¬ 
pending on Cyg X-1 number Ncygx-i numerical modelling 
results are presented in the two figures. Note that we cal¬ 
culate number of BH-l-Psr systems as the next ratio - 
1500 • NsH+Psr/Npsr, where 1500 is rounded number of 
known radio pulsars, Nbh+Pst is the number of BH-|-Psr 
binaries and Npsr is the total number of radio pulsars ap¬ 
peared during population synthesis. 

There were no Cyg X-1 type systems in all of the mod¬ 
els with type B wind including Wb model, so all of these 
scenarios must be declined right away because of we observe 
the object and do not regard it as a statistical ejection. 

In the Figure 1 models for the mass loss scenarios A 
and C at a very wide changing of all of the unknown pa¬ 
rameters are shown. Since we know only one Cyg X-1 type 
X-ray source (namely Cyg X-1) and purpose of this work is 
to display that under any suppositional parameters BH-|-Psr 
systems have to be observable, we have presented in Figure 
1 and Figure 2 only such models in which number of Cyg 
X-1 type systems is no more than 4 and number of BH-|-Psr 
binaries is no more than 3. ’’Forbidden zone” in these hg- 
ures is highlighted area where there are no any models with 
conceivable appropriate parameters. 

Black squares (mod.l. Figure 1) denote the model cal¬ 
culated with kbh = 0.45, minimal black hole mass Mmin = 
7 Mq, mass loss type A and uniform initial mass ratio distri¬ 
bution {aq = 0). Characteristic kick velocity vo for neutron 
stars and parameter for black holes take on a values 0, 
180, 360, 720 km s“^. We change these velocities for NS and 
BH both to give a demonstration of it’s influence on depen¬ 
dence between BH-fPsr and Cyg X-1 quantities. Number of 
BH+Psr systems if uo = 0 and Vq^ = 0 is much more than 3, 
so this point is not presented. Models in which characteristic 
kick velocity for neutron stars and parameter for black 
holes take on a values 720, 360, 180 km s“^ are shown in 
order of increasing of Cyg X-1 type systems number (higher 
kick reduces of number of binaries including relativistic com¬ 
panion) . 

Grey circles (mod.2, Figure 1) depict models with stel¬ 
lar wind type A, uniform initial mass ratio distribution 
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Figure 1. Dependence between BH+Psr binaries and Cyg X-1 type systems quantities. Black squares (mod.l) denote the model 
calculated with = 0.45, minimal black hole mass Mmin = mass loss type A and uniform initial mass ratio distribution 

{^OLq = 0). Characteristic kick velocity vq for neutron stars and parameter for black holes take on a values 0, 180, 360, 720 km s~^. 
Number of BH+Psr systems if uq =0 and Vq^ = 0 is much more than 3, so this point is not presented. Models in which characteristic kick 
velocity for neutron stars and parameter Vq^ for black holes take on a values 720, 360, 180 km s“^ are shown correspondingly in order 
of increasing of Cyg X-1 type systems number. Grey circles (mod.2) depict models with stellar wind type A, uniform initial mass ratio 
distribution {aq = 0), characteristic kick velocity vq = 180 km s“^ for neutron stars and parameter Uq^ = 180 km s~^ for black holes. 
During calculations quantity k^h has been varied within the limits 0.1 and 1 every 0.1, black hole minimal mass takes on a values 2.5, 3.0, 
4.0, ..., 9.0, 10.0 Mq. Circles in the Figure 1 evidently group in some vertical lines. Each line corresponds to one value of k^hi minimal 
black hole masses change along lines. Cyg X-1 systems usually are products of evolution of very massive binaries, so difference between 
their number in case of Mmin = IOMq and in case of Mmin = 2.5Mq is negligible. Otherwise number of binary radio pulsars with 
black holes has strong dependence on minimal black hole mass - the more Mmin the less BH+Psr systems number. So the bottom-up 
sequence order of points in lines which depict models with various is: Mmin = 10? 9,..., 4, 3, 2.5Mq. Parameter kf^^ influences on Cyg 
X-1 and BH+Psr numbers both, the more k^h, the more number of binaries including relativistic companion. Note that in many cases 
quantities of Cyg X-1 or BH+Psr systems are higher than limits in Figure 1. That is why in this figure presented only that models in 
which kbh ^ 0.6 (for mod.2) and not all points in vertical lines. So in mod.2 vertical lines correspond (in order of decreasing NcygX — i) 
to the next k^p values: 0.6, 0.5, 0.4, 0.3. Models with k^h = 0.1, 0.2 has no Cyg X-1 type binaries and are merged with many other 
models without such systems in this figure. Open squares (mod.3) describe models with stellar wind type C and quadratic initial mass 
ratio distribution (aq = 2). Quantity minimal black hole mass Mmin-, characteristic kick velocity for neutron stars vq and parameter 
for black holes in mod.3, mod.4 and mod.5 take on the same values as in mod.2. Also in these models as in mod.2 points group into 
some vertical lines corresponding their k^h-, minimal BH mass change along lines. Maximum k^h for the presented models (mod.3) is 
0.7. Black triangles (mod.4) mark models with stellar wind type A and quadratic initial mass ratio distribution (aq = 2). Maximum k^h 
for the presented models is 0.6. Black diamond formations (mod.5) designate models with stellar wind type C and uniform initial mass 
ratio distribution (aq = 0). Maximum k^h for the presented models is 0.6. ’’Forbidden zone” is highlighted area where there are no any 
models with conceivable appropriate parameters of stellar evolution. 


(aq = 0), characteristic kick velocity vq = 180 km s“^ 
for neutron stars and parameter Vq^ = 180 km s“^ for 
black holes. During calculations quantity khp bas been var¬ 
ied within the limits 0.1 and 1 every 0.1, black hole minimal 
mass takes on a values 2.5, 3.0, 4.0, ..., 9.0, 10.0 Mq. Cir¬ 
cles in the Figure 1 evidently group in some vertical lines. 
Each line corresponds to one value of kbh , minimal black hole 
masses change along lines. Cyg X-1 type systems usually 
are products of evolution of very massive binaries, so differ¬ 


ence between their number in case of Mmin = lOM© and 
in case of Mmin = 2.5M0 is negligible. Otherwise number 
of binary radio pulsars with black holes has strong depen¬ 
dence on minimal black hole mass - the more Mmin the less 
BH+Psr systems number. So the bottom-up sequence or¬ 
der of points in lines which depict models with various kbh 
is: Mmin = 10, 9,..., 4, 3, 2.5M0. Parameter kbh influences 
on Cyg X-1 and BH+Psr numbers both, the more kbh, the 
more number of binaries including relativistic companion. 
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Note that in many cases quantities of Cyg X-1 or BH+Psr 
systems are higher than limits in Figure 1. That is why in 
this figure presented only that models in which kbh ^ 0.6 
and not all points in vertical lines. So for mod.2 vertical 
lines correspond (in order of decreasing Ncygx-i) to the 
next khh values; 0.6, 0.5, 0.4, 0.3. Models with k^h ~ 0.1, 0.2 
has no Cyg X-1 type binaries and are merged with many 
other models without such systems in the Figure 1. 

Open squares (mod.3, Figure 1) describe models with 
stellar wind type C and quadratic initial mass ratio distribu¬ 
tion {aq = 2). Quantity kbh, minimal black hole mass Mmin, 
characteristic kick velocity for neutron stars vq and param¬ 
eter for black holes take on the same values as in mod.2. 
Also as in mod.2 points group into some vertical lines cor¬ 
responding their kbh, minimal BH mass change along lines. 
Maximum kbh for the presented models is 0.7. 

Black triangles (mod.4, Figure 1) mark models with 
stellar wind type A and quadratic initial mass ratio dis¬ 
tribution (aq = 2). Quantity kbh, minimal black hole mass 
Mmin, characteristic kick velocity for neutron stars vo and 
parameter for black holes take on the same values as 
in mod.2. Also as in mod.2 points group into some verti¬ 
cal lines corresponding their kbh, minimal BH mass change 
along lines. Maximum kbh for the presented models is 0.6. 

Black diamond formations (mod.5, Figure 1) designate 
models with stellar wind type C and uniform initial mass ra¬ 
tio distribution(a 5 = 0). Quantity kbh, minimal black hole 
mass Mmin, characteristic kick velocity for neutron stars vo 
and parameter for black holes take on the same values 
as in mod.2. Also as in mod.2 points group into some verti¬ 
cal lines corresponding their kbh, minimal BH mass change 
along lines. Maximum kbh for the presented models is 0.6. 

As one can see from the Figure 1 in the value area con¬ 
taining at least one black hole candidate being in Cyg X-1 
type system all of the models are complying with the condi¬ 
tion 1500- (Nbh+Pst) /Npsr ^ 0.4. In zone which has limits 
1 ^ Ncygx-i ^ 3 lower limit of number of binary radio 
pulsars with black holes is between « 0.4 and « 1.75 in all 
models. 

In the Figure 2 two kinds of Wc models are shown (Wb 
model is declined because it has no Cyg X-1 type systems). 

Upright crosses (mod.6, Figure 2) describe model with 
uniform initial mass ratio distribution (uq = 0). During cal¬ 
culations characteristic kick velocity vq for neutron stars and 
parameter Vq^ for black holes take on a values 0, 180, 360, 
720, 1000 km s“^. If uo = 0 and Vg^ = 0 quantity of Cyg X-1 
systems and BH-(-Psr binaries are more than limits of the 
Figure 2 and appropriate point is not presented. So upright 
crosses depict mod.6 and crossbucks depict mod.7 in which 
characteristic kick velocity for neutron stars and parameter 
Vq^ for black holes take on a values 180, 360, 720, 1000 km 
s“^ are shown correspondingly in order of decreasing of Cyg 
X-1 type systems number. 

Crossbucks (mod.7) denote model with quadratic ini¬ 
tial mass ratio distribution (oq = 2). During calculations 
characteristic kick velocity vg for neutron stars and param¬ 
eter Vg^ for black holes take on a values 0, 180, 360, 720, 
1000 km s“^. If Do = 0 and Vg^ = 0 quantity of Cyg X-1 
systems and BH-fPsr binaries are more than limits of the 
Figure 2 and appropriate point is not presented. So upright 
crosses depict mod.6 and crossbucks depict mod.7 in which 
characteristic kick velocity for neutron stars and parameter 



Figure 2. Dependence between BH+Psr binaries and Cyg X-1 
type systems quantities (Wc model). Upright crosses (mod.6) de¬ 
scribe model with uniform initial mass ratio distribution (aq = 0). 
Crossbucks (mod. 7) denote model with quadratic initial mass ra¬ 
tio distribution (aq = 2). During calculations characteristic kick 
velocity vg for neutron stars and parameter Vg^ for black holes 
take on a values 0, 180, 360, 720, 1000 km s"*^ in both cases. If 
uq = 0 and Vg^ = 0 quantity of Cyg X-1 systems and BH+Psr 
binaries are more than limits of the Figure 2 and appropriate 
point is not presented. So upright crosses depict mod.6 and cross¬ 
bucks depict mod.7 in which characteristic kick velocity for neu¬ 
tron stars and parameter Vg^ for black holes take on a values 180, 
360, 720, 1000 km s”*^ are shown correspondingly in order of de¬ 
creasing of Cyg X-1 type systems number. ’’Forbidden zone” is 
highlighted area where there are no any models with conceivable 
appropriate parameters of stellar evolution. 

Vg^ for black holes take on a values 180, 360, 720, 1000 km 
s“^ are shown correspondingly in order of decreasing of Cyg 
X-1 type systems number. 

It is evidently from the Figure 2 that in the value area 
containing at least one black hole candidate being in Cyg 
X-1 type system all of the models are complying with a 
condition 1500 • (NpH+Par)/Npar 0.4. In zone which has 
limits 1 ^ Ncygx-i ^ 3 lower limit of number of binary 
radio pulsars with black holes is between ~ 0.4 and ss 1.75 
in all models. 


5 CONCLUSIONS 

iLipunov et al.l (ll994^ made a conclusion that observ¬ 
able by modern techniques BH+Psr binaries have to be in 
Galaxy. Despite new theories concerning evolution of mas¬ 
sive stars appearance we confirm this conclusion about a 
possibility to discover binary radio pulsars with black holes. 
We suggest that expected value of pulsars paired with black 
holes comparative abundance may be found within the lim¬ 
its 0.4 < . 1500 < 2. We also confirm the conclu¬ 

sion made by TLinunov et all ^1994^ about high eccentricity 
of such systems (Figure 3) and about sufficient number of 
system close enough to observe (Figure 4). Distribution of 
eccentricities shows evident peak at e ~ 1 that is the con- 
sequenc e of mass loss and kick dur ing second supernova ex¬ 
plosion iKornilov fc Lipunovlll98^ . 

Our calculations are important for estimation of 
merging of BH+BH and BH+NS systems. Results of 
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Figure 3. Eccentricity distribution of BH+Psr binaries in Wc 
model with quadratic initial mass ratio distribution {aq = 2), 
characteristic kick velocity for neutron stars vq and parameter 
for black holes take on a value 180 km s“^. 

Number of systems, % 



Figure 4. Orbital period distribution of BH+Psr binaries in Wc 
model with quadratic initial mass ratio distribution {aq = 2), 
characteristic kick velocity for neutron stars vq and parameter 
Ug^ for black holes take on a value 180 km s“^. 

iLipunov fc Panchenkol ll20nff> shows that merging rate can 
increase by « 5 — 7 times with respec t to previous compu¬ 
tations made bv iLipunov et all New estimations of 

merging rate of relativistic systems will be carried out in one 
of our future work. 

In closing we emphasize again that results of this work 
does not depend on optimal evolutionary scenario param¬ 
eters: radio pulsars paired with black holes have to be in 
Galaxy and might be discovered within the next few years. 
We accentuate that relative number of BH+Psr systems to 
the total number of pulsars have been calculated in this 
work and most of the pulsars do not experience recycling 
effect. Radio pulsars in BH+Psr binaries originated in our 
calculations have no differences with more than 90% of radio 
pulsars which have already been observed. 
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Table 1. Possible evolutionary track which produces BH+Psr. 


Stage of 
primary 
star (1) 

Stage of 
secondary 
star (2) 

Mi/Mq 

M 2 /M 0 

a/RQ 

Time 

10®yr 

CE-stage 

MS 

MS 

69.5 

5.1 

700 

0.0 

- 

SG 

MS 

66.4 

5.0 

730 

3.3 

- 

Rov 

MS 

63.5 

5.0 

760 

3.6 

J- 

WR 

SN 

Rov 

37.9 

5.1 

13 

3.6 

— 

BH 

Rov 

16.4 

4.9 

410 

3.7 

- 

BH 

WR 

SN 

16.4 

4.9 

14 

4.9 

- 

BH 

Psr 

16.4 

1.34 

21 

26.3 

- 


Table 2. 

Possible evolutionary track which produces Psr+BH. 

Stage of 
primary 
star (1) 

Stage of 
secondary 
star (2) 

Mi/Mq 

M2jMQ 

a/RQ 

Time CE-stage 
10®yr 

MS 

MS 

23.9 

11.5 

450 

0.0 

SG 

MS 

22.4 

11.4 

480 

6.3 

Rov 

MS 

21.0 

11.4 

500 

6.9 

WR 

MS 

8.5 

23.9 

690 

6.9 

SN 






Psr 

MS 

1.34 

23.9 

790 

7.1 

Ej 

SG 

1.34 

22.3 

830 

13.0 


SN 





Psr 

BH 

1.34 

9.0 

80 

14.0 
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APPENDIX: POSSIBLE EVOLUTIONARY 
TRACKS WHICH PRODUCE BH+PSR 
BINARIES 

We have presented possible evolutionary track which 
produces BH+Psr binary in Table 0 and possible evolu¬ 


tionary track which produces Psr+BH system in Table 
Note that they do not depict all possibilities and evolu¬ 
tion of concrete binary strongly depends on evolutionary 
scenario. In both cases we used stellar wind type A, param¬ 
eter khh = 0.43. 

Marks in the tables Q and |21 depict the next stages: 
MS - main sequence stage, SG - super-giant stage, Rov - 
Roche lobe overflow stage, WR - Wolf-Rayet star stage, BH 
- black hole, Psr - radio pulsar, Ej - ejecting neutron star 
which does not appear itself as radio pulsar due to free- 
free absorbtion of radio emission in component’s stellar wind 
(detailed nome nclature of neu tron stars may be found in 
book written bv lLiDun*^ EsS)): SN - supernova explosion. 

First and second columns in each table contain informa¬ 
tion about primary (more massive) and secondary compan¬ 
ion current stage correspondingly, third and fourth - about 
their masses, fifth contains value of major semi-axis of the 
system, in sixth we have presented time elapsed from the 
moment of binary system birth. Seventh column answer the 
question - is the system in common envelope stage or not? 

Tracks which produce neutron star before black hole 
can lead to so-called ’’recycled” radio pulsars forming. In 
general, we calculate evolution of rotation of a neutron star 
and can include these pulsars into consideration. Neverthe¬ 
less, in this case results of our work would be dependent on 
many extra assumptions: distribution of neutron stars on 
their magnetic held, magnetic held dissipation time, inhu- 
ence of accretion on neutron star’s magnetic held, etc. We 
prefer not to use any additional hypothesis since part of such 
tracks is negligible (see Introduction). 










